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Numerous experiments on the detection of atmospheric waves in the 
frequency range from acoustic to planetary at meteor heights have revealed 
that important wave sources are meteorological processes in the troposphere 
(cyclones, atmospheric fronts, jet streams, etc.)- A dynamical theory 
based on the works of OBUKHOV (1949) and MONIN (1958) include describing 
the adaptation of meteorological fields to the geostropic equilibrium 
state. According to this theory, wave motions appear as a result of 
constant competition between the maladjustment of the wind and pressure 
fields due to non-linear effects and the tendency of the atmosphere to 
establish a quasi-geostrophic equilibrium of these fields. 

Barotropic Atmospheric Model 

To demonstrate the above, we consider an approximation of the 
atmosphere to be a two-dimensional liquid film located in the Coriolis 
force field. The equations for this model are obtained by averaging over 
height the equations of motion and continuity for a three-dimensional 
atmosphere (OBUKHOV, 1949). Then, passing over to dimensionless variables, 
we choose values corresponding to large-scale synoptic motions (OBUKHOV, 
1949) b^ way of the scales of length L and velocity W. The natural scale 
(2 a) ) corresponding to OBUKHOV's (1949) characteristic adaptation time 
will serve as the temporal scale. In the dimensionless variables we have 
equations : 
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where 2u) is a Coriolis parameter; U, V ? are the wind velocity components; C, 
= x/L; n Z = y/L; U = u/W; V = v/W; n = c x/2w WL; 

t = 2co z t ; x = In ( p/p Q ) ; 

p is pressure, p is a standard pressure near the earth. The dimensionless 
parameters of Rossby-Kibel R = W/2to L and p = 2w L/c have been introduced 
into (1). To analyze the system (1), we use a small Rossby number for 
motions of the synoptic scale and apply the method of asymptotic series and 
multiple time scale successive approximations (BLUMEN, 1972; JEFFREY and 
KAWAHARA, 1982), according to which we introduce, apart from the "fast" 
dimensionless time t = 2to t, the "slow" dimensionless time T = R x = tw/L 
and change 
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The zero-order system in Ro 
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was investigated by OBUKHOV (1949). OBUKHOV showed that a complete 
solution of the nonlinear system (2) is the superposition of the non 
stationary waves and the stationary geostrophic components. 


In the first-order approximation in Ro the system (1) takes the form: 
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In the analysis of wave generation within the flux under the influence 
of nonlinear terms we consider the atmosphere in a quasi-adjusted state 
when the process of adaptation either has been completed and the wave 
component has already been scattered in space or the wave motion has not 
occurred at all due to the initial adjustment of the meteorological fields. 
Therefore, in what follows, we shall consider the zero-order terms in (3) 
to be purely geostrophic and independent of the "fast" time x. Proceeding 
in the usual way we obtain from (3) equations for the vorticity 
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By excluding from the third equation of (3) and the second equation of 
(4), we obtain 
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In the right-hand side of (5) there is a function T, in the left-hand one 
there are functions x and T. Therefore, to eliminate divergence in (5), 
according to the requirements of the method of multiple time scales, we 
assume that the right- and left-hand sides of (5) separately equal zero. 

By equating to zero the right-hand side of (5) we obtain the equation of 

slow quasi-geostrophic evolution of a potential vortex (BLUMEN, 1972). 
Then^ from the condition of the left-hand side of (5) being zero, we obtain 
11^ = 0. Thus, the first order terms in Ro in the ejyansion (5) are 
purely wave terms. Making use of the condition £1^ - (5 11^ = 0, it is 

possible to derive wave equations with a nonlinear forcing term on the 

right-hand side. For instance, the equation for has the form: 
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The expression on the right-hand side of (6) describes the rate of wave 
generation by the background quasi-geostrophic motions. This process 
exists owing to the nonlinear transfer of energy between motions of 
different types: slow quasi-geostrophic and fast waves. 


Influence Functions of Tropospheric Meteorological Sources 

We can obtain an equation for the spectral density of the vertical 
velocity W (GAVRILOV, 1985) for an isothermal three-dimensional atmosphere: 


+ n 2 W = F U, w, k , k ); £ 2 = c/H; k 2 = kx 2 + k 2 ; 

d £-2 w x y y 

2 ,..2 2 W , 2 2, . 2 u 2 , 2 u 2 lr 2 ... 

n = (N - to ) / ( oj - t ) k H + (joH/C - 1/4; 

s 


F = 


w 




its) exp ( -C/2 ) ; 


(7) 


where N is the Brunt frequency; W, Q, and S are the Fourier transforms of 
the vertical components of wave velocity W and the source terms (4) q and 
s, respectively; w, k , k are the frequency and horizontal projections of 
the wave vector, respectively . Analysis of (7) shows that the contribution 
of the source Q increases in the region of high frequencies. For the 
acoustic frequencies ou >> 2co the main contribution to the generation of 
acoustic waves is made by Q, the expression for which coincides with that 
for the source of the sound generation by turbulence obtained by LIGHTHILL 
(1952). For the low frequency IGW the contribution of the addendum 
containing s is essential whereas the values of Q and S are smaller and of 
the same order. 

Using the condition of radiation W ~ exp (iqu) at u -» °° and being 
interested in a solution for large heights only, one can find the solution 
of (7) by the method proposed in KAMKE, (1976) and GAVRILOV (1985): 
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Here the wave field of the vertical velocity W is the Fourier integral: 
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In (9) we proceed from integration over spectral variables Q, k , k to 
integration over the natural variables t, x, y. With this in view, w^ use 
the theorem of the Four ier-transf orm of the convolution of functions and 
find the influence functions of the sources q and S: 
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The expressions for W through and G^ have the form: 
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where the symbol 6 denotes the convolution of functions in the variables 
x, y, t. It can be seen from (11) that the wave amplitudes increase 
exponentially with height. Therefore, of interest is the behavior of G^ 
and G^ at u + To calculate the integral 

exp[i(n£ - k x - k y)] dk dk , (12) 

K x y x y 

contained in (10) we use an asymptotic expansion by the stationary phase 
method (FELSEN and MARCUVITZ, 1973) at C -*■ 00 . We obtain, for example, for 
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where y = R/(2HN); R = + y 2 ; C* = 4R 2 u> 2 + N 2 H 2 . 

z c, 

The expression for is derived in a similar way but is omitted here 

because of its bulkiness. The first addendum in (13) describes the high 
frequency mode of IGW with the frequencies to 'v N, and the second addendum 
is the lower frequency waves with periods of the order of hours. 

The described procedure of obtaining the Green's function was used in 

the calculation of the wave field of the vertical velocity W for a typical 
meteorological formation, such as the cyclonic vortex. The mathematical 
model of the latter is given by the expression 

p = p Q - 6p*exp(- t 2 /t 2 - x 2 /a 2 - y 2 /b 2 ); (14) 

where the cyclonic pit depth p changes with height 
p ( u ) = p n exp[-(z - z g ) 2 /h 2 ] 

In (14) p is the standard pressure at a given height; a, b are the scales 
of the vortex length along the corresponding axes, x is the vortex 
characteristic lifetime; z is the vortex maximum height; is the vertical 
halfwidth of the source. S In the calculations the values which correspond 
to mesoscale synoptical vortices of the moderate type were used: a = 50 

km, b = 100 km, h = 2 km, 6p m = 20 m bar, Z g = 8 km, x = 3 hrs. Fig. 1 

shows the directional pattern of radiation in the plane x, y for the 

vertical velocity W. It can be seen that the maximum radiation intensity 
takes place along the semimajor axis of the vortex. The ratio of the 
radiation intensities along the main semiaxes increases in proportion to 
the square of the a/b ratio. Fig. 2 shows a horizontal cross-section along 
the x axis of the wave field of the vertical velocity W. It can be seen 
that the vertical velocity values reach tens of cm s , and the decrease of 
the amplitude is due to the source attenuation with time. The temporal 
variation of W at the point x = 1,500 km, y = 0.3, H = 100 km is presented 
in Fig. 3. It can be seen that oscillations at a given point occur with a 
period T ~ 80 min and are modulated by a wave with the period T ~ 16 hrs. 
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